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Abstract. In this paper we describe some analytical properties of the 
^i/m ^g^^ proposed by Sersic (1968) to categorize the photometric pro- 
files of elliptical galaxies. In particular, we present the full asymptotic 
expansion for the dimensionless scale factor h{ra) that is introduced when 
referring the profile to the standard effective radius. Surprisingly, our 
asymptotic analysis turns out to be useful even for values of m as low 
as unity, thus providing a unified analytical tool for observational and 
theoretical investigations based on the P}/^ law for the entire range of 
interesting photometric profiles, from spiral to elliptical galaxies. The 
systematic asymptotic analysis provided here also allows us to clarify the 
value and the limitations of the power law often used in the past 
as a natural representation of the surface brightness profiles of elliptical 
galaxies. 



1. Introduction 

After its introduction as a generalization of the i?^/^ law (de Vaucouleurs 1948), 
the i?^/™ law (Sersic 1968) has been widely used in observational and theoretical 
investigations (see, e.g., Ciotti & Bertin 1999, and references therein, hereafter 
CB). According to such law, the surface brightness profile is given by 

I{R) = he-^"^"'^ , m>0, (1) 

where r] = R/Re, and 6 is a dimensionless constant such that Re is the effective 
(half-luminosity) radius. The projected luminosity inside R is given by 

L{R) = 2vr I{R')R'dR' = I^Rl^^{2m, br,'/^), (2) 

where 7(0, x) = e~^t'^~^dt is the (left) incomplete gamma function. The total 
luminosity is then given by 

L = IoRf-^Ti2m), (3) 
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where r(a) = 7(0, 00) is the complete gamma function. It follows that b{m) is 
the solution of the following equation: 

7(2™,6) = nM. (4) 

2. Asymptotic expansion 

Unfortunately, eq. (4) cannot be solved in explicit, closed form, and so it is 
usually solved numerically. This is inconvenient for a number of observational 
and theoretical applications. At least three interpolation formulae for b{m) 
have been given in the literature, namely 6 ~ 1.9992m — 0.3271 by Capaccioli 
(1989), 6 - 2m - 0.324 by Ciotti (1991), and 6 - 2m - 1/3 + 0.009876/m by 
Prugnicl &: Simien (1997). These expressions provide an accurate fit in the range 
0.5 < m < 10; curiously, their leading term is linear in m, with a slope very 
close to 2. In CB it is proved that this behavior results from a general property 
of the gamma function. In particular, it turns out that the first terms of the 
asymptotic expansion of 6(m) as implicitly given by eq. (4) are 

'^^^ 4 + 4oL + 255tL? + 

Equation (5) now clearly explains the value of the interpolation formulae found 

earlier. Note that 4/405 = 0.009876 Of course, the asymptotic analysis 

provided in CB would allow us to give explicitly any higher order term if so 
desired. In CB it is shown that this expansion, even when truncated to the first 
four terms as in eq. (5), performs much better than the interpolation formulae 
cited above, even for m values as low as unity, with relative errors smaller than 
10-^ 

The use of this simple formula is thus recommended both in theoretical and 
observational investigations based on the Sersic law. As a simple application, in 
CB the leading term of the asymptotic expansion of the total luminosity, of the 
central potential, and of the surface brightness profile associated with the R^/™' 
law are obtained in explicit, analytical form. In particular, the connection of eq. 
(1) with the simple power law R'^, often used in the past to fit the photometric 
profiles of elliptical galaxies, is brought out explicitly. 
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